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A characterization of qubit quantum channels is introduced. In analogy to what happens in the 
context of Bosonic channels we exploit the possibility of representing the states of the system in 
terms of characteristic function. The latter are functions of non-commuting variables (Grassmann 
variables) and are defined in terms of generalized displacement operators. In this context we intro- 
duce the set of Gaussian channels and show that they share similar properties with the corresponding 
Bosonic counterpart. 
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In quantum mechanics the transition from the initial 
state to the final state of a system is described in terms of 
quantum channels [H . At a mathematical level these are 
linear maps operating on the set of bounded operators 
of the system, which preserve the trace and (if any) the 
positivity of the operators on which they act. Finally 
in order to represent a "physical" transformation, i.e., 
a transformation that could be implemented in a lab, a 
quantum channel must also possess the property of com- 
plete positivity (i.e., the positivity of any initial joint op- 
erator acting on the system plus an external ancilla need 
to be preserved by the action of the map). An impres- 
sive effort has been devoted in the last decades to study 
the properties of quantum channels. Indeed they play 
a fundamental role in many different branch of physics, 
specifically in all those sectors where one is interested in 
studying the decoherence and noise effects. 

In the context of quantum information theory empha- 
sis is put on characterizing the properties of quantum 
channels in terms of their information capacities 0, Q . 
These figures of merit are the quantum counterpart of the 
Shannon capacity of a classical communication line 
which "measure" the performances of the map in con- 
veying classical or quantum information. Even though 
impressive achievement has been obtained in this field in 
the recent years, several open questions are still under 
investigation — we refer the reader to Q and references 
therein for details. 

The majority of the results obtained so far relate to 
two specific classes of channels, namely the qubit chan- 
nels and the Bosonic Gaussian channels. The former are 
completely positive trace-preserving (OPT) transforma- 
tions which act on the state of a single two-level quantum 
system (qubit). Due to the small size of the Hilbert space 
a simpleparametrization of these channels has been ob- 
tained [a, 0] while some additive issues jl, la, lOl and 
several classical and quantum capacities [s Uglllll lia. [l3| 
have successfully been solved (see also Ref. Q for a re- 
view) . Bosonic Gaussian channels [3, [IB| > on the con- 
trary, are a specific subclass of OPT maps acting on a 
continuous variable system that preserve certain symme- 
tries. These channels include a variety of physical trans- 
formations that are of fundamental interest in optics, in- 
cluding thermalization, loss and squeezing. As in the 



qubit channel case, additivity issues |16l . Il7| and capaci- 
ties 18, id, HO, HH, nil have been successfully solved for 
Bosonic Gaussian channels. Furthermore, they allow for 
a compact parametrization [13, [Hi [23l [23 in terms of 
the characteristic function formalism 251 2a , [2?! . 

In this paper we establish a parallelism among the 
qubit channels and the Bosonic Gaussian channels by 
introducing for the former a characteristic function rep- 
resentation. To do so we adapt the formalism introduced 
by Cahill and Glauber in Ref. [1^ for representing the 
density operators of Fermions to the case of two-level 
systems. In this context the channels are represented in 
terms of Green functions. Interestingly enough this al- 
lows us to define a set of Gaussian channels for qubit that 
share analogous properties with their continuous variable 
counterpart. 

The paper is organized as follows. In Sec. U we briefly 
review the characteristic function formalism for Bosonic 
(and Fermionic) systems. In Sec. [Ill we introduce the 
displacement operator and characteristic function for a 
qubit. To do so we introduce Grassmann variables and 
we use them to generalize the definition of coherent states 
for finite dimensional systems. We then present a Green 
function representation for qubit channels (Sec. IIIip and 
define the set of qubit Gaussian channels (Sec. IIVP dis- 
cussing their degradability properties. Conclusion and 
final remarks are presented in Sec. IVl The paper in- 
cludes also a couple of technical Appendixes: namely, 
in Appendix VK\ we review some properties of Grassmann 
calculus, while in Appendix [Bl we present a brief excursus 
on quantum channel degradability. 



I. CHARACTERISTIC FUNCTION FOR 
BOSONS 

In quantum optics a complete description of the state 
of a Bosonic mode characterized by annihilation and cre- 
ator operators a and a^, can be obtained in terms of its 
coherent states These vectors possess various ap- 

pealing properties. Specifically, they minimize the un- 
certainty relations of any couple of conjugate quadratures 
and they are eigenvectors of a. Most importantly, coher- 
ent states form an over-complete continuous set of vec- 
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tors parametrized by a single complex variable fi. This 
allows us to expand any other state of the system as a 
superposition of the |/i)s with coefficients which define 
quasi-probability density functions. Exploiting this and 
the fact that the coherent states can be obtained by ap- 
plying the displacement operator D{fi) = exp[/ia^ — /i*a] 
to the vacuum, one can also use the latter as an over- 
complete operator basis |14', 25, 27]. In particular, given 
any bounded operator of the system (e.g., a density 
matrix p), we can write 



where d?^ — (iRc(/i)rfIm(/i) and where 
x(At) = Tr[ei?(Ai)] . 



(1) 



(2) 



Equation ([2]) defines the characteristic function of the 
operator Q. This is a complex function of the variables 
/i and fi* which provides us with a faithful description 
of the original operator thanks to the "orthogonality" 
relation 



TT[Di^,)D{-,y)]=6('\^,~,y) 



(3) 



with (5*^^^(/i — ly) being the Dirac delta in the complex 
plane. To represent a density matrix p, the function ([2]) 
needs to possess certain properties [Tj, |23 including be- 
ing continuously differentiable in /i and /i* and verifying 
x(0) = 1. Within the characteristic function description, 
Gaussian states are defined as those p whose x(m) is a 
Gaussian function of the complex parameter /i (exam- 
ples arc thermal, coherent, and squeezed states). 

Consider now the action of a linear super-operator TV 
which transforms 6 into 9' = A/'(6). Equation ([T]) allows 
us to represent this mapping in terms of a linear trans- 
formation of x(m)- Indeed the characteristic function of 
the output operator Q' is 

= / d^,^xMG{,^,fi) , (4) 



with 



G(i.,p) = TT[D{-i^)MH[D{fi)y 
= Tr[AA(z?(-^))z?(/i)]. 



(5) 
(6) 



In these expressions Mh is the dual of TV which describes 
the channel in the Heisenberg picture and which is de- 
fined by the identity 



Tr[eiAA^(e2)] = Tr[AA(ei)e2] , 



(7) 



for all Oi.2 (see, for instance, [2l[). We call Eq. (O the 
Green function of TV: according to the previous defini- 
tions it provides us with a complete characterization of 
the channel. 



A special subset of OPT map s for Bosonic systems is 
the set of Gaussian channels [l^ . These are characterized 
by Green functions ^ of the form 



Giv^p) — (J*-^-* (i^ — wp — ui/i*) 

xexp -i(M*, -M)r 



(8) 



with r being a real symmetric positive matrix (i.e., 
covariance matrix) and v and w being complex num- 
bers — rigorously speaking, Eq. ([8]) defines one-mode 
Bosonic Gaussian channels. As can be directly verified 
from Eq. Q such maps have the peculiar property of 
transforming input Gaussian states into output Gaussian 
states. An interesting fact about these channels is that, 
except for the additive classical noise channel [2l|, l23l |. 
they admit a physical representation j2ll | in terms of a 
single mode environment originally prepared in a Gaus- 
sian state. Indeed, the exceptional role of the additive 
classical noise channel corresponds to the fact that any 
one-mode Bosonic Gaussian channel can be represented 
as a unitary coupling with a single-mode environment 
plus an additive classical noise. Within such represen- 
tation (without additive classical noise) one can show 
that the Bosonic Gaussian channels ([8]) are either anti- 
degradable or weakly degradable [l^]. Moreover, in the 
case in which the single-mode representation is of Stine- 
spring form (that is, if the environment state is pure) 
the channel is then anti-degradable or degradable in the 
sense of Ref. (9j| (for the sake of completeness explicit 
definitions of these properties are given in Appendix |B]) . 



A. Characteristic function for Fermions 

The characteristic function formalism presented in the 
previous section can be generalized to describe Fermionic 
systems too [1^. The main difference in this case is 
related to the fact that now the complex variables /i 
and p* are replaced by a couple of conjugate Grass- 
mann variables ^ and £,* [1^ whose properties are re- 
viewed in Appendix [X] This is intrinsically related with 
the fact that the annihilation and creation operators of 
a Fermion obey anti-commutation rules instead of com- 
mutation rules [sO]. We will not review the analysis of 
Ref. ^28*1 since in the next section, when discussing the 
qubit case, we will rederive most of the results obtained 
in the Fermionic case. 



II. REPRESENTATION OF A QUBIT 

Various proposals for defining a (discrete) phase space 
for finite dimensional systems have been discussed so far 
by introducing generalized position and momentum oper- 
ators (see, for instance, Ref. f3l| and references therein). 
Here we will not follow this line: instead we invoke the 
analogies between a qubit and a single Fermionic mode to 
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adapt the results of Ref. [28[ . A similar approach was de- 
veloped in Ref. [33] to solve non-Markovian master equa- 
tions of a two-level atom interacting with an external 
field. 

The starting point of our analysis is to observe that the 
lowering and raising operators of the qubit [i.e., (7+ = 
|1)(0| and (T_ = satisfy anti-commutation rules 

similar to those of a Fermionic mode, i.e., 

{r7_,a+} = |0)(0| + |1)(1| = 1, 

= {a+,(7+} = 0. (9) 

Identifying the qubit state |0) with the Fermionic vacuum 
we can therefore treat <7+ and tT_ as Fermionic creation 
and annihilation operators, respectively. Following (28j 
we introduce then a couple of conjugate Grassmann vari- 
ables ^ and ^* (see Appendix [X| and impose standard 
anti-correlation with the annihilation and creator opera- 
tors of the system, i.e., 

u,^±} = {r,^±} = o. (10) 

It is worth noticing that this implies that the projectors 
0)(0| ~ (T_(T+ and |1}(1| = cr+cr_ as well as the Pauli 
matrix Oz = |0)(0| — |1)(1| commute with f and 
In the following we will also require that 

eij) = (-i)^IjK 

rij) = (-i)^ij)r, (11) 

for j = 0, 1. This is not strictly necessary but it is con- 
sistent with Eq. and allows us to simplify the calcu- 
lations. For instance, given any collection of qubit oper- 
ators 8i, 82, • • • , On+i and the Grassmann numbers ^1, 
^2, • • • , n we can use Eq. (jlip to verify that the following 
relation applies 

Tt[ei(ie2^2---enCnOn+i] (12) 

(an analogous expression holds also when replacing all, 
or part of, the with their complex conjugates — more 
details about the trace can be found in Appendix I A 1|) . 

The above definitions give us the possibility of oper- 
ating with "hybrid" mathematical objects obtained by 
multiplying Grassmann variables and qubit operators. In 
this context we find it useful to define a generalized ad- 
joint operation for these hybrid operators by arbitrarily 
imposing the conditions 

(616626 •••0nene„+i)t 

= 0i+ie:et... 6*eUi*0l , (13) 

with 6 and Qi as in Eq. (fT2ll . 

A. Qubit characteristic function 

Qubit displacement operators can now be defined in 
analogy with [2^ as 

D{0 = exp(a+e-6^-) (14) 



where in the second line we used Eq. (|A6p . As in the 
Bosonic case they satisfy the identity -D^(6 = 
Moreover the application of -D(6 to the vacuum orig- 
inates eigenvectors of the annihilation operator of the 
system (cr_). These are the coherent states of our qubit, 
i.e., 

I0 = ^(6|o) = (i-^)|o)-6i), (15) 

whose norm is unity. These vectors are eigenvectors of 
(T_ in Grassmann sense (i-e^ their eigenvalues are Grass- 
mann variables; see Ref. [28| for details). 

What is interesting for us is the fact that D{^) can be 
used to define a characteristic function for the operators 
of the system as in Eq. i.e., 

x(6 ^ neDiO] . (16) 

In particular, consider a 6 which is characterized by 
the matrix 

when expressed in the computational basis {|0), |1)}. In 
this case using the anti-commutation rules of Eq. (fTO|) 
and the identity (fT^ we get 

X(6 = Tr[e] + {9oo - Oii)^ + SoiC - ^106- (18) 

It is worth noticing that with respect to the analysis of 
Ref. ^] the characteristic functions analyzed here con- 
tain an extra term which is linear in ^ and Indeed 
in the Fermionic case analyzed by Cahill and Glauber 
the only allowed physical states are classical mixtures 
of |0)(0| and |1)(1| (this follows from the requirement 
of invariance under 27r rotation with respect to an ar- 
bitrary axis). Consequently the off-diagonal terms asso- 
ciated with ^01 SLiid 9iQ do not need to be considered. 
When analyzing qubit systems, instead, quantum super- 
positions among |0) and |1) are allowed and we need to 
include also the linear contributions. 

As in the Bosonic case, Eq. p6)) can be inverted. In 
this case, however, Eq. ([T]) is replaced by 

e = y d'e xiO E{-0 , (19) 

with E{^) ^ D{^) defined by 

E{0 = a, - 6e/2 + <J+C - Ct- . (20) 

The easiest way to verify this is by direct substitution of 
Eqs. (HH) and into Eq. HH) and by employing the 
integration rules (jA7|l . 
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1. Density operators 

To represent a density operator 
_ ( p 7 



P ■ 



7* 1 



(21) 



the characteristic function needs to satisfy certain phys- 
ical requirements. First of all, the Hermitianity of p and 
the normalization condition Tr[p] = 1 imply, respectively. 



x(o) - 1, 



(22) 
(23) 



where complex conjugation is defined as in Eq. (|A4|) [to 
verify this simply use Eq. (fT8|) with O = p]. The posi- 
tivity of p imposes, instead, the following inequality to 
hold 



1 



(24) 



This follows from the positivity condition I7P ^p{l~p) 
and by the identity 



7 - j d^£. x(.o r , 

d'i xiO + 1/2 



P 



Using similar arguments one can verify that Eqs. ([22]) - 
(|24| are also sufficient conditions for x(C) being a char- 
acteristic function of a density operator p. 



III. GREEN FUNCTION REPRESENTATION 
OF A QUBIT CHANNEL 

Let us now consider the effect of a qubit quantum chan- 
nel M acting on the operator 8 of the system. As in the 
Bosonic case we would like to derive its Green function 
representation ([6]) . To do so we first evaluate the charac- 
teristic function x'iO associated with AQB with A and 
B being arbitrary qubit operators. This is 

= I d\ TT[AxiOEi-OBDiO] , (25) 



where we used Eq. with xiO being the characteristic 
function of (from now on C and ^ should be considered 
entries of the same Grassmann set). Our goal is to find 
a function G(C, £,) which gives 



d\xiO G(C,6 



(26) 



for all xiO- Notice that if ^ were a commuting variable 
(e.g., a complex variable) the problem could be solved by 



simply moving xiO out of the trace operation of Eq. ([25]) 
yielding G{(,0 = Ti[AE{-C)BD{£_)]. In the case under 
consideration, however, the situation is complicated by 
the fact that for moving out of trace the variables ^ or ^* 
we need to insert cr^s as in Eq. (jl2p . Taking into account 
this fact, the solution becomes 



C?(C,0 = TT[Aa,Di-C)BD{0] , 



(27) 



as can be easily verified by direct integration of the 
Eqs. and (|26p for the most general characteristic 
function 

The Green function (|26p associated with a CPT map TV 
can then be obtained by using an operator sum represen- 
tation [l|, Q of such channel and exploiting the linearity 
of the trace. Indeed, writing A/'(6) = J2k ^'^k^Ml with 
{Mk}k being Kraus operators of TV, we get 



k 

= TrI 



Using Eq. (|A13p this can also be written as 
G(C,6 = Trla^Di-CWnlDiO 



(28) 



(29) 



with Afn being the Heisenberg representation of the map 
TV defined in Eq. ([7]). Equation (^9)1 shows that, as in 
the Bosonic case, a complete description of the channel is 
obtained by applying the dual map to the displacement 
operator — see Eq. Exploiting the normalization 

condition J^k ^l-^k = 1 we note that for ^ = the 
above expression yields 



G(C,0) = Tr[a,i^(~C)] =CC 



(30) 



which corresponds to the Grassmann delta function 
(5'^^(C) defined in Eq. (jA9p . in agreement with the 
requirement of channel being trace preserving — see 
Eqs. and (pS )) . 

Finally, let TVi and TV2 be two different qubit channels 
with Green functions Gi(C,0 G2{C,£,)j respectively. 
From the definition ([^5]) we then find that the Green func- 
tion Gi2(C, of the composite map TV2 oTVi in which we 
first operate with TVi and then with TV2 , can be expressed 
in terms of the following Grassmann convolution integral 

GuiCO = / d'e Gi(C,C') G2(e',0, (31) 
with C, ^, and ^' Grassmann numbers. 

A. Examples and canonical forms 

As a particular case of Green function consider the 
identity map X which leaves all operators invariant, i.e., 
X(8) ~ 8. According to our definition we get 

Gico = TA'^M-OD{o] = ic-oic-n , (32) 
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which, as expected, corresponds to the delta 5^^^(C — 
of Eq. (|A9[) . More generally from Ref. Q we know that 
the most generic qubit quantum channel TV implements 
the following transformation. 



M{p) = M 



1 + r-a 



l + {t + Tr) ■ 3 



(33) 



where t — {t\^ti^t-i) is a real vector, a = {ui, 172, cts} 
is a vector containing the Pauli matrices, r is the Bloch 
vector describing the input state, and T is a real 3x3 
matrix. In Ref. 6] it is shown that T can be reduced, 
via changes of basis in (i-e., via proper rotations of 
the input and output states), to the diagonal (canonical) 
form T — diag{Xi, X2, A3), with the real coefficients 
^^1,2,3 and ii,2,3 that need to satisfy certain conditions 
0] to guarantee the complete positivity of the map. In the 
Green function language such canonical form corresponds 
to have 



(34) 



^(2)(^_^i±^^_^i^^*,exp 



+(A3 - AiA2)er + ^i-^cce - ^^^ccr • 



IV. GAUSSIAN CHANNELS FOR QUBITS 

In analogy with the Bosonic case, in this section we 
introduce the definition of qubit Gaussian channels. We 
start noticing that in order to define these channels it 
does not make sense to focus on maps which transform 
Gaussian characteristic functions into Gaussian charac- 
teristic functions. Indeed, thanks to Eq. (jA6p . all charac- 
teristic functions of a qubit can be written in a Gaussian 
form Therefore following Eq. ([5]) we say that a 

qubit map is Gaussian if its Green function has the form 

G(C, = ^^'^ (C - - bC) exp[-cre] , (35) 

with a and b complex and c real [s^ l numbers, respec- 
tively, and with the exponential defined as in Eq. (jA6p . 
A trivial example is provided by the identity map I whose 
Green function ([5^ is of the form ([55)1 for 6 = c = and 
a = 1. 

Generic mixtures of Gaussian channels do not neces- 
sarily have the form ([55)1 . Therefore the set of Gaus- 
sian channels is not convex. However, it has semi-group 
structure with respect to the channel composition rule o. 
Indeed, given two Gaussian channels TVi and ^2 charac- 
terized by parameters (ai,5i,ci) and (02,62,22), respec- 
tively, from Eq. ([5T)) it is easy to verify that the Green 
function of Af2 o A/i is again of the form ()35p with 

a aia2 + 6162 , 
b = 0162 + hal , 

|2 



Cl(|a2r - I62I ) +C2 



(36) 



Both the semi-group property and the non-convexity 
property hold also in the Bosonic case. 



Canonical form for Gaussian channels 



From Eq. ([34)l it is easy to verify that within the 
parametrization 0] we can get Gaussian maps ()35p 
by choosing 



A3 

h 



A1A2 , 

t2^0 



(37) 
(38) 



This in fact yields Gaussian Green functions with a — 
(A2 + Ai)/2, 6 = (A2 - Ai)/2 and c = ^3/2. We can then 
use [7] to show that the corresponding transformation is 
CPT if and only if the following inequalities hold. 



\Xk\ ^ 1 forfeit, 2; 

\t,\ ^ V(1-AD(1-A^). 



(39) 



This allows us to parametrize the whole set of Gaussian 
channels in terms of three real parameters only. First 
of all, as in Refs. 0, we can use a trigonometric 
parametrization to express Ai_2 in terms of the angles 9, 
(f) in [0, 2ti[ as follows 



Ai = cos( 



A2 = cos( 



(40) 



Then we can parametrize ^3 by introducing the positive 
quantity q G [0, 1] to write 



H - (2g - 1 



cos(26') - cos{2(j)) 



(41) 



Replacing all this into Eq. ([51)1 yields the following canon- 
ical form for the Green function of a qubit Gaussian chan- 
nel, i.e., 



G(C,0 =^^^^ (C-Ccos6icos 
X exp 



^* sin 9 sin c 



_ cos(2g)-cos(20) ^^^^. 



(42) 



We will see that the maps of this form have the peculiar 
property that they can always be described in terms of 
a unitary interaction of the form ()Bip with a single (not 
necessarily pure) qubit environment. For this reason we 
call them "qubit-qubit" channels. It is worth stressing 
that once again a similar property holds for the Bosonic 
case: there (almost) all the one-mode Bosonic Gaussian 
maps are in fact describable in terms of a single mode 
environment [l^, [2l| . 



B. Qubit-qubit maps: Pure environment case 

An important subclass of the qubit-qubit channels of 
Eq. is obtained for q = 1 and 9 and generic, i.e., 

G(C, = ^^^^ (C - C cos 61 cos -f C sin 9 sin 0) 
cos(2(?) — cos(2(/ 



X exp 



(43) 
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According to Eq. pQ]) this corresponds to having Its] = 
- Xl){l - Xl). As shown in Ref. 0] any CPT map 
which can be described in terms of an interaction with 
a single qubit environment originally prepared in a pure 
state can be expressed in this form by proper unitary 
rotation of the input and the output state. This implies 
that the maps (|43p admit a Stinespring dilation (|Bip with 
a two-dimensional (qubit) environment E. Without loss 
of generality, we can assume an initial state of the envi- 
ronment of the form pE = |0)£;(0|. Following Ref. 
one can then choose the unitary coupling U to have the 
following block structure 



U = 



[Ai] [a^Aoa^] 



with 

Ao 



cos 9 
cos ( 



Ai 




sin( 



sm( 




(44) 



(45) 



being a Kraus set for the channel [the matrix (|44| is 
expressed in the basis {|00), |10), |01), with \jk) = 
\j)(E)\k)E ioT J, k = 0,1]- 

The complementary channel TV" [s^ [s^ can now 
be computed as in Eq. (jB2p . Since it represents a 
qubit channel — it connects two two-dimensional Hilbert 
spaces (the input Hilbert space with the environmental 
one) — we can use Eq. (|28p to evaluate its Green function 
obtaining 



G{C,£,) = 5(2) (C-Ccoseisin^ + f sin6icos( 
"cos(26i) -hcos(20) 



X exp 



(46) 



It is still of the (pure-environment qubit-qubit) Gaussian 
form (|43p and can be expressed in terms of the original 
Green function G{(, £,) of A/" by simply shifting (p by — 7r/2 
and by changing sign to 9, i.e., 



-e 



(47) 



In Ref. 11^ it has been shown that qubit-qubit 
channels with pure environment are degradable for 
cos(26')/ cos(20) ^ 0, and anti-degradable otherwise. 
Here we will rederive this same result in the Green func- 
tion formalism as a consequence of the Gaussianity of 
these maps, pointing out an interesting parallelism with 
their Bosonic counterpart. 

In analogy with [lo, HH we look for the intermediate 
map T that should connect M with J\f, in the class of 
qubit-qubit channels (with pure environment). Rewrit- 
ing the degradability condition (jBSp in terms of the com- 
positions rules ((5T|) we can then recast the problem as 
follows 



(48) 



where GxiCO is the Green function (|43p of the map T 
characterized by the parameters 9^ and By using 
Eq. (136|) we find that, for cos(26')/ cos(20) ^ 0, 9^,(l3x 
do exist such that Eq. is satisfied. Specifically such 
parameters are defined by the relations. 



COs(20a;) 

cos(2(/)2,) 



cos(26') - cos(20) -f 2 008(26*) cos(20) 

cos(26l) + cos(20) 
cos(26i) - cos(20) - 2cos(26l) cos(20) 



cos(26l) + cos(20) 



(49) 



The case cos(20)/ cos(20) ^ can be treated analo- 
gously to show that the corresponding channels are anti- 
degradable. In fact, in the Green function formalism the 
anti-degradability condition (|B4p becomes 



(50) 



where Gx{C,S.) is the Green function of the connecting 
map r. We find that for cos(26l)/ cos(20) 0, Eq. ^ 
is satisfied by choosing Gx{C,0 ™ the subclass of qubit- 
qubit channels with pure environment - i.e., Eq. - 
with 9x and (j)x determined by the expressions (HS)) after 
replacing {9, (f) with (—0, (j) — 7r/2). 

More directly this result can be established by using 
the correspondence ((47|) and the fact that the comple- 
mentary channels of degradable maps are anti-degradable 
— see Appendix IbI Consider, in fact, a (pure environ- 
ment) qubit-qubit channel N with 008(26*)/ cos(20) ^ 0. 
According to Eq. (|47)) we know that its complemen- 
tary M is still a (pure environment) qubit-qubit chan- 
nel characterized by the parameters [9' — {~9,(p — 
7r/2). Now it is easy to verify that cos(20')/ cos(20') — 
-cos(26l)/cos(2(/)) ^ 0. Therefore from Eqs. (gH]) and 
we can conclude that Af is degradable while TV is 
anti-degradable. 

Note that, in the special case 008(26*) — cos(20) — 0, 
both the degradability relations are satisfied. Therefore 
in this case the qubit-qubit channels with pure environ- 
ment are both degradable and anti-degradable, with null 
quantum capacity. 



C. Qubit-qubit maps: Mixed environment case 



Now let us consider the Gaussian channels (|42)) for 
q ^ 1. They can be represented in terms of a physical 
representation (jBip with U as in Eq. (|44p and with E 
being a single qubit environment initially prepared in the 
mixed state, 



PE = q\0)E{0\ + il-q)\l)E{l\ 



(51) 



To verify this, we observe that with the above prescrip- 
tions Eq. (|Bip gives 

(52) 



J^{p) - TrE[U{p^[q\0)E{0\ + {l 
= qN^ip) + {I - q)Mi{p) , 
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with Mo = Tte[U{p^ \0)e{0\)W] bei ng the (pure en- 
vironment) qubit-qubit channel of Sec. IIVBI associated 
with the operator U and with Mi{p) = a-xM'o{(Jxpcrx)(^x- 
From the properties of ax it follows that a Kraus set for 
Afi is given by the matrices (|45p by exchanging 9 and (/). 
Consequently the Green function of this channel is given 
byG(C,C)|e^0 withG'(C,<^) asinEq. (gS]). Using this fact 
and the linear dependence of Eq. (|28p with respect to TV 
we can now evaluate the Green function of the map (152]) 
as follows 

G(C,0 = g(5*^^ (C-?cos6lcos(/) + rsin6'sin(/)) 
cos(26') — cos(2(? 



X exp 



+ (l-q) S^^'> (C-^cos^cosei + rsin^sine*) 
"cos(2(?!)) - cos(26') 



X exp 



(53) 



Equation ([55]) can finally be casted into the form 
thanks to the identity 



-X cc* 



l + {2q-l)x^C 

J2q~l)xie 



(54) 



which holds for all x complex — see Eq. (jA6p . The above 
is an example of a convex combination of Gaussian chan- 
nels (i.e., Afo and A/i) which is still Gaussian. 

A natural question is then whether or not the weakly 
complementary channel (jB2p associated with Eq. ([5^ 
is also Gaussian. To see this we first use the linearity 
of trace to express the complementary A/" as a convex 
combination of the weakly complementaries of A/q and 
TVi, i.e., J\f = q Afo + (l - q) j(fi. Then we invoke the 
linearity of Eq. (^5)) and use Eq. (|^ to write 



"cos(26') -Fcos(2(^) 



X exp 



+ (l-g) (^(2) (C + ?sin(/)cos6l-rcos(/)sin6') 
cos(2^) + cos(26') 



X exp 



(55) 



This is of the form (|42|) only for q = 0, 1. Therefore, in 
general, the weakly complementaries of qubit-qubit maps 
with mixed environment are not Gaussian even though 
they can be expressed as a convex combination of Gaus- 
sian channels (i.e., Afo and A/i). This can be pushed a 
little further by observing that for generic choices of 0, 
4> and q, the weakly complementaries (1551) are not even 
unitarily equivalent to a qubit Gaussian channel [37| . 



1. Weak-degradability properties 

Let us analyze the weak-degradability properties of the 
qubit-qubit channels with mixed environment. 



As in Sec. IIVBI we prove that the maps Af of Eq. (142]) 
are weakly degradable for cos(20)/ cos(20) ^ 0. In this 
regime in fact one can easily check that Eq. (|48|) can 
still be solved with Gx{C,£,) of the form (|55p replacing 9 
and (j) with —9x and (jfx + 7r/2 where Ox,4'x satisfy the 
relations (f49l) . 

Proving anti-degradability for cos(20)/ cos(20) ^ is 
not simple because, in general, Af is not in a Gaussian 
form — see Eq. ([55]) . However, in this case we show 
that these channels cannot be used to transf er q uantum 
information since their quantum capacity Q [38| is null. 
To see this we notice that for cos(20)/ cos(2(/)) < 0, TV is 
a mixture l|52|) of two channels (i.e., Afo and TVi) which 
are both anti-degradable and have hence null quantum 
capacity, i.e., Q{Afo) = Q{Afi) =0 — see Appendix IB] 
Under these conditions it is easy to verify that also TV 
must have a null Q. Indeed let us consider a new CPT 
map, 

N'{p) = q AAo(p) ® |0)b(0| + (1 - g) A/'i(p) ® |1)b(1| , 

where B is an ancillary system. We can now verify that 
the TV is isomorphic to f o TV' with — TrB[...] ® 

0)_b(0| being a CPT map which replaces all states of B 
with a fix given output |0)b- Expressing Q in terms of 
the output coherent information |39| of the channel and 
using the quantum data processing inequality |l| we can 
verify that Q{Af) ^ Q(TV'). Besides, by using the basic 
properties of von Neumann entropy [l| we can express the 
coherent information of TV' as J(TV', p) = qJ{Afo, p) + {l — 
q)J{Afi,p). Putting all this together we get 

QiAf') = lim maxJ([A/'']«^,p)/Af 

< gQ(A/'o) + (l-g)Q(A/'i) -0, (56) 

and hence Q{Af) — 0. 



V. CONCLUSIONS 

In this work we introduce a characteristic function for- 
malism for the qubit channels in terms of generalized dis- 
placement operators and Grassmann variables, inspired 
by a parallelism among these maps and the Bosonic 
Gaussian channels. 

We then present a Green function representation of 
the quantum evolution that allows us to define the set 
of qubit Gaussian maps. In this context, we find that 
all the Gaussian channels are qubit-qubit, i.e., they can 
always be described in terms of a unitary interaction of 
a qubit system with a single (not necessarily pure) qubit 
environment. Similarly, it is known that in the Bosonic 
case (almost) all the one-mode Bosonic Gaussian maps 
are describable in terms of a single mode environment. 

This formalism turns out to be elegant and powerful 
and, in particular, it can be used to study the weak- 
degradability properties of the qubit-qubit maps, for 
both pure and mixed qubit environments, in terms of 
Green functions. 
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On one hand, in the case of pure environment, the 
qubit-qubit maps are either degradable (i.e., additive co- 
herent information) or anti-degradable (i.e., Q=0). Be- 
sides, the complementary maps are still qubit-qubit chan- 
nels and so Gaussian. It is interesting to note that an 
equivalent property holds for one-mode Bosonic Gaus- 
sian channels. On the other hand, in the case of mixed 
environment, we show that the qubit-qubit maps are ei- 
ther weakly degradable or they cannot be used to transfer 
quantum information (i.e., Q=0). However, in this case 
the weakly complementary maps do not belong to the set 
of qubit-qubit channels and are not Gaussian. 

It is important to stress that this Green function for- 
malism shows clearly that the qubit Gaussian maps share 
analogous properties with their continuous variable coun- 
terpart, i.e., the Bosonic Gaussian channels. 

Finally, we remark that the characteristic function ap- 
proach, introduced in this paper for qubit systems, can be 
generalized to d-levcl quantum systems (qudit) in terms 
of generalized Grassmann variables [401 . 
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APPENDIX A: GRASSMANN VARIABLES 

A Grassmann variable <^ spans over a set of objects (the 
Grassmann numbers) ^i, ^2, which anti-commute. 
Indeed, given any and elements of the set, they 
satisfy the relation 



C»6 



(Al) 



while obeying ordinary commutation relations with re- 
spect to the multiplication by a complex number. In 
particular Eq. (jAip implies that a Grassmann variable 
is 2-nilpotent, i.e., = (note that is trivially in- 
cluded in the Grassmann variable set). At a mathemati- 
cal level, the above conditions can be rigorously formal- 
ized by saying that Grassmann numbers are the gener- 
ators of an algebra over the complex field which obey 
anti-commutation relations. 

Complex conjugation of ^ can be defined by introduc- 
ing an extra Grassmann variable ^* whose elements ^l, 
^2, ••■ obey the same relation (|Aip and anti-commute 
with all the ^^s, i.e.. 



CO = -oc 



(A2) 
(A3) 



To identify with the complex conjugate of we finally 
require the relations 



to be satisfied for any x complex number or product of 
thea,6,--- and^i*,^^,---. 

Given the above properties it follows that the most 
general function /{£,,£,*) is linear both in ^ and i^*, i.e., 



(A5) 



with A, Bi,2, and C independent from ^ and In 
particular, the exponentials become 



exp(i3ie + + ceo = E 



n=0 

= 1 + + B2C + ere + Bi^B2C/2 + B2CBi^/2 . 

(A6) 

This expression can be used to verify that (apart from 
a global multiplicative term) any function (|A5|) can be 
written as an exponential. 

Integration over ^ and ^ can be defined by introduc- 
ing the "differential" and ■ These are assumed to 
obey the same anti-commutation relations obeyed by the 
variables ^ and including Eqs. ([Aljl . (fTO|) . and (fTTjl . 
The integrals are then defined according to the Bcrczin 
rules 



dc r = 1 



(A7) 



Joint integration with respect to ^ and ^ is finally de- 
fined by identifying the double differential d^^ as follows. 



d^£,=dCd^ = -d^dC 



(A8) 



In this context one can identify an analogous of the Dirac 
delta function S^'^^ (/i — i') in the complex plane. Such 
Grassmann delta is defined as 

S(^\^~0 ^ Jd\exp[KiC~C)~{^~0^^*] 

= (e-c) (r-c) , (A9) 

with e, C, and k Grassmann variables. Indeed, from 
Eq. (jA7|) and from Eq. (jASp we have 



(AlO) 



5('He-c)/(e,r)-/(c,c), 



for aU fiOO- Notice that the delta function (|A9|) com- 
mutes with any Grassmann numbers and satisfies the re- 
lation s<-^\i - c) = s^^^c -0 = -s'^^Hc - O- 

A useful property is the following. Given the function 
/{£,,£,*) one can define its even and odd parts, i.e., 

(AU) 

According to Eq. (|A5|1 they are of the form /+{£,,£,*) = 
A + and /-(CD = Bi ^ + B2 T, respectively. 

Now given g{^,S.*) another function we can write 



(A4) 



I d'e/±(e,r)5T(e,r 



= 
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and thus 



+ fd'u-i^,n9A^.n- (A12) 



1. More about trace 

Equation shows that the cychcity of the trace 
needs to be modified when involving Grassmann terms. 
If we need to move only qubit operators, then the stan- 
dard rule applies, i.e., 

iv[ei6 • • ■ e„5„e„+i] = Tr[e„+ieiei---e„en] 
= Tr[a---e„Cne„+iei] . 

(A13) 

On the contrary, if we move also Grassmann variables, 
by exploiting the anti-commutation rules of the ^^s, we 
get 

Tr[eiae26---e„6^e„+i] (au) 
= (-i)"-iTr[c„e„+ieiei---e„] 
= (-i)"-iTr[e26 • • • e„e.©n+iei6] . 

Finally in conjunction with Eq. (|13p. Eq. (TT^ gives 

(Tr[eiae26 • • • e„Cne„+i])* (ai5) 



APPENDIX B: WE AK-DEGRAD ABILITY VS. 
ANTI-DEGRAD ABILITY 

It is a well known (see, e.g., [4l|, [i^) that any CPT 
map TV can be described by a unitary coupling between 
the system S with an external ancillary system E (de- 
scribing the environment) prepared in some fixed pure 
state. This follows from the Stinespring dilation 43] of 
the map which is unique up to a partial isometry. More 
generally, one can describe TV as a coupling with an en- 
vironment prepared in some mixed state pE-, i.e.. 



where Tr £;[...] is the partial trace over the environment 
E and U is a unitary operator in the composite Hilbert 
space TLs®'He- As proposed in Ref. [2l| we call Eq. (jBl[) 
a "physical representation" of TV to distinguish it from 
the Stinespring dilation, and to stress its connection with 
the physical picture of the noisy evolution represented by 
TV. Moreover, Eq. (|B1[) motivates the following defini- 
tion [i^, [2lj . For any physical representation in Eq. (|B1[) 
of the quantum channel TV we define its weakly comple- 
mentary as the map TV which takes the input state p into 
the state of the environment E after the interaction with 
S, i.e., 



N{p)^i:is[u{p®pe)u^] 



(B2) 



The transformation (|B2p is CPT and describes a quan- 
tum channel connecting systems S and E. It is a gener- 
alization of the complementary channel TVcom defined in 
Refs. [1, [sH, [3^. If some channel T docs exist such that 



{roN){p)=M{p) , 



(B3) 



for all density matrices p, then M is called weakly degrad- 
able and J\f anti-degradable. Similarly if 



(ToAA)(p)=AA(p), 



(B4) 



Afip)^TYE[U{p(g)pE)U^ , 



(Bl) 



for some channel T and all density matrices p, then Af 
is anti-degradable while A/" is weakly degradable (see [23, 
[2TI'). In Ref. [S] the channel TV is called degradable if 
one considers the environment in a pure state. Clearly 
any degradable channel [qI is weakly degradable but the 
opposite is not necessarily true. 

Degradability and anti-degradability have been proved 
useful to analyze the quantum capacity [ssl of the chan- 
nel. On one hand, one can verify that anti-degradable 
channels (where this property is defined irrespectively 
from the purity of pE associated with the physical repre- 
sentation) cannot be used to convey quantum messages 
in reliable fashion — i.e., their quantum capacity Q nul- 
lifies [H, [2^, [21I . On the other hand, instead degradable 
channels [9| allows for a single letter formula expression 
for Q — i.e., the maximum of their output coherent in- 
formation is additive. 
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